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Abstract 

Quantum decoherence, the evolution of pure states into mixed states, may be a feature 
of quantum-gravity models. In most cases, such models lead to fewer neutrinos of all active 
flavours being detected in a long baseline experiment as compared to three-flavour standard 
neutrino oscillations. We discuss the potential of the CNGS and J-PARC beams in con- 
straining models of quantum-gravity induced decoherence using neutrino oscillations as a 
probe. We use as much as possible model-independent parameterizations, even though they 
are motivated by specific microscopic models, for fits to the expected experimental data 
which yield bounds on quantum-gravity decoherence parameters. 
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I. INTRODUCTION 



If microscopic black holes, or other defects forming space-time foam, exist in the vacuum 
state of quantum gravity (QG) , this state, in our view, will constitute an "environ- 
ment" which will be characterised by some entanglement entropy, due to its interaction with 
low-energy matter. Based on earlier work by Maldacena jsj] on anti-de-Sitter space black 
holes Hawking [4| has recently claimed the absence of information loss in quantum gravity. 
He argued that in a superposition of topologies, the non-unitary contributions associated 
with the non-trivial topology decay exponentially with time, leaving only contributions to 
the path integral from the unitary topologically-trivial configurations. However we do not 
believe that this settles the issue of induced decoherence in quantum gravity since firstly the 
concept of anti-de Sitter space times cts db rc gulator" has been used and secondly a Euclidean 
formulation of the path integral for quantum gravity has been adopted. Both these features 
may not be shared by an eventual true theory of quantum gravity. Furthermore, entangle- 
ment entropy is still present for an outside observer (la] and space-time foam need not consist 
only of black holes; other defects, e.g. point-like D-branes, might be present. A stochastic 
fluctuation of populations of point-like D-particles 6|, for instance, which are known to obey 
infinite statistics [7] due to their infinite internal stringy states, could constitute such a deco- 
herening environment for matter propagation. There are also other quantum-gravity issues 
that are not completely understood. In particular, it is possible that the entire concept of 
local effective Lagrangians breaks down in such situations. 

The matter system in such a case behaves as an open quantum mechanical system, 
exhibiting decoherence, which has in principle detectable experimental signatures. In the 
context of a phenomeno logical parametrization of quantum-gravity induced decoherence the 
first tests along these lines have been proposed in [8( . A more microscopic consideration was 
given in [9|, where the proposed parametrization of decoherent effects of quantum gravity 
was forced to obey the Lindblad 10, lUl] formalism of open systems, employing completely 
positive dynamical semigroup maps. This latter phenomenology, however, may not be a 
true feature of a quantum theory of gravity. 

The decoherent approach to quantum gravity, entailing entanglement entropy, has been 



followed by some of the authors 



12. 



13| in many phenomenological tests or microscopic 



models of space-time foam [6[], within the framework of non-critical string theory; the latter 
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may be a viable (non-equilibrium) theory of space-time foam [141 ]. based on an identification 
of time with the Liouville mode. The latter is viewed as a dynamical local renormalizat ion- 
group scale on the world-sheet of a non-conformal string. The non-conformality of the 
string is the result of its interaction with backgrounds which are out of equilibrium, such 
as those provided by twinkling microscopic black holes in the foam. The entropy in this 
case can be identified with the world-sheet conformal anomaly of a cx-model describing the 

n 

propagation of a matter string in this fluctuating background [L4J . Although within critical 
string theory, arguments have been given that entanglement entropy can characterise the 



number of microstates of Anti-de-Sitter black holes [IJjj], we do not find these to be entirely 
conclusive, and moreover an extension of such counting to microscopic dynamical black holes, 
that characterise a space time foam, is far from being understood. For instance, the process 
of formation and annihilation of microscopic black holes or other singular fluctuations in 
space time, including defects, may not be described by critical string theory methods. 

In view of the above issues, it is evident that the debate concerning space-time foam 
remains open. The thermal aspects of an evaporating black hole are suggestive that the 
environment due to quantum-gravity is a sort of "thermal" heat bath. This has been pursued 
by some authors, notably in ref. [if]]. Another proposal, the D-particle foam model {(J, 
considers the gravitational fluctuations that could yield a foamy structure of space-time 
to be D-particles (point-like stringy defects) interacting with closed strings. There are no 
thermal aspects but there is still the formation of horizons and entanglement entropy within 
a fluctuating metric framework. 

In general, for phenomenological purposes, the important feature of such situations is the 
fact that gravitational environments, arising from space-time foam or some other, possibly 
semi-classical feature of QG, can still be described by non-unitary evolutions of a density 
matrix p. Such equations have the form 

d t p = Aip + A 2 p (1) 

where 

AiP = \ [ P , H] 

and H is the hamiltonian with a stochastic element in a classical metric. Such effects may 



arise from back-reaction of matter within a quantum theory of gravity [8|, 113] which deco 



heres the gravitational state to give a stochastic ensemble description. Furthermore within 



models of D-particle foam arguments in favour of a stochastic metric have been given 12]. 
The Liouvillian term A 2 p gives rise to a non-unitary evolution. A common approach to 
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11] but this is not 



A 2 p, is to parametrise the Liouvillian in a so called Lindblad form 
based on microscopic physics. We note at this point that any non-linear evolutions that may 
characterise a full theory of QG (see e.g. a manifestation in Liouville strings 181]). can be 
ignored to a first approximation appropriate for the accuracy of contemporary experimental 
probes of QG. Generically space-time foam and the back-reaction of matter on the gravita- 
tional metric may be modelled as a randomly fluctuating environment; formalisms for open 
quantum mechanical systems propagating in such random media can thus be applied and 
lead to concrete experimental predictions. The approach to these questions have to be phe- 
nomenological to some degree since QG is not sufficiently developed at a non-perturbative 
level. 



One of the most sensitive pro 



neutrinos 



19 



20 



21 



22 



23 



:>es of such stochastic quantum-gravity phenomena are 



24 



25l . l26l . l27j . and in particular high-energy ones [28, 



291 ] . It is the point of this article to study decoherence induced by non-linear space-time 



foam fluctuations as a subdominant effect in neutrino oscillations at CNGS and J-PARC 
beams after giving an overview of the framework of decoherence phenomena in neutrino 
experiments. 



A linear decoherence simplified model, of Lindblad type [IOJ has been used for the fit, 
following earlier work in 30]. The model of 3(| involved a diagonal decoherence matrix, and 
in the fit of (which we will refer to as I), decoherence was assumed to be dominant only 
in the antineutrino sector, in order to fit the LSND results 31| pointing towards significant 

—* oscillations, but suppressed oscillations in the particle sector. In this way in I 
a fit was made to a three generation model with the LSND "anomalous" result, without 
introducing a sterile neutrino. The strong CPT violation in the decoherence sector, allowed 
for an equality of neutrino mass differences between the two sectors, so as to be in agreement 
with atmospheric and solar neutrino data. However, the reader should bear in mind that 
the LSND result has not been finally confirmed by MiniBoone 32|. Nevertheless, this does 
not detract from the motivation to use a decoherence-based formalism to discuss quantum 
gravity effects for neutrinos. 

The particular choice of I, which yielded the best fit to all available neutrino data in- 
volved mixed energy dependence for the (antineutrino-sector) decoherence coefficients, some 
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13|, 



of which were proportional to the neutrino energies E, while the rest had been inversely 
proportional to it, oc 1/E. In I, the coefficients proportional to 1/E were interpreted as 
describing ordinary matter effects, whilst the ones proportional to E were assumed to cor- 
respond to genuine quantum-gravity effects, whose increase with the energy of the probe 
was consistent with the fact that the higher the (anti)neutrino energy the larger the back 
reaction effect on the quantum space time, and hence the larger the decoherence. The strong 
difference assumed in I between the decoherence coefficients of the particle and antipartic 
sectors, although not incompatible with a breakdown of CPT at a fundamental level 
appears at first sight somewhat curious; in fact it is unlike any other case of decoherence for 
other sensitive particle probes, like neutral mesons, examined in the past js,!^]. There, the 
oscillations between particle and antiparticle sectors, necessitate a common decoherence en- 
vironment between mesons and antimesons. If one accepts the Universality of gravity, then, 
the best fit of I seems incompatible with this property. Moreover, there are two more prob- 
lematic points of the fit in I, which were already discussed in that reference. The first point 
concerns the complete positivity of the model. In [30|] the diagonal form of the decoherence 
matrix, used in I, was taken ad hoc, without explicit mention of the necessary conditions to 
guarantee complete positivity, as required by the Lindblad approach. Furthermore, the par- 
ticular choice of the best fit of I, did not lead to positive definite probabilities for the entire 
regime of the parameter space of the model, although the probabilities were positive definite 
for the portion of the parameter space appropriate for the various neutrino experiments used 
for the fit. Specifically, it was found in I, that with the particular choice of the decoherence 
parameters in the (antineutrino sector), one obtains positive-definite transition probabilities 
for energies E > 0(1 MeV). The second, and more important point, is that the best fit of 
I is good for all the neutrino experiments available at the time, but unfortunately it could 
not reproduce the spectral distortion observed by the KamLand experiment [34j , whose first 
results came out simultaneously with the results of I. In (which we will refer to as II) 
the above weakness was rectified by requiring general conditions among the coefficients that 
guarantee complete positivity in the entire parameter space for the three-generation simpli- 



fied Lindblad linear model of decoherence of 



301 ] . used in I. In fact it was shown in II that 



it was possible to choose the coefficients in a way to give a consistent and excellent fit to all 
available neutrino data including the spectral distortion seen by KamLand, and the LSND 
results 31[ for the transition probabilities in the antineutrino sector. 
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It has been argued that quantum decoherence could be an alternative description 



of neutrino f 
Kamiokande 



avour transitions. Fits to data by different experiments such as Super- 



35| and KamLand 34( have been performed and these clearly disfavour a 



decoherence explanation for neutrino oscillations. However, quantum decoherence may still 
be a marginal effect in addition to neutrino oscillations and could give rise to damping fac- 
tors in the transition probabilities reducing the number of active neutrinos being detected 
in a long baseline experiment compared to what is expected from the standard neutrino 
oscillation scenario. 

Our article is organized as follows. In section [Til we describe the physics of the "quantum- 
gravitational analogues" of the MSW effect and of foam models endowed with stochastic 
fluctuations of the space-time metric background. In section IIII[ we discuss decoherence 
signatures in the neutrino oscillations and review the existing bounds on these parameters 
using the available neutrino data, including those from KamLand 34(, indicating spectral 
distortions. Then, in section IIVI we present the damping signatures and the associated 
fitting functions, which might be due to either the "quantum-gravitational analogue" of the 
MSW effect or the stochastic fluctuations of the space-time metric background. We are 
careful to consider various stochastic models of foam, which lead to different dampin g sig na- 
tures, depending on the details of the underlying characteristic distribution functions [36| . In 
section [V] we estimate the sensitivity of CNGS and J- PARC experiments to the parameters 
of quantum-gravitational decoherence entering the set of the above-mentioned damping sig- 
natures. Finally, in section IVIl we compare the sensitivities estimated by means of bounds 
obtained from data on atmospheric, solar and KamLand neutrino oscillations, as well as 
the neutral kaon system, and discuss the possible relevance of our results in guiding the 
construction of models of (the still elusive theory of) Quantum Gravity. 



II. THEORETICAL MODELS FOR QUANTUM- GRAVITY DECOHERENCE 
AND NEUTRINOS 

The picture we envisage is the following: there are several parallel three-brane worlds, 
one of which represents our observable Universe, embedded in a higher- dimensional bulk, in 
which only gravitational (closed) string states propagate. On the brane world there are only 
open string states propagating, representing ordinary matter, with their ends attached on 
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the hypersurface. Of course, there are also closed string states, either propagating along the 
longitudinal brane directions, or crossing the brane boundary from the bulk. As discussed 



in 



371 ] consistent supersymmetric models of D-particle foam can be constructed, in which 



the bulk space between, say, two parallel brane worlds is populated by point-like D-particle 
defects. Motion of either D-particles or branes, as required by the need to have cosmological 
backgrounds for a brane observer, causes supersymmetry breaking in both the brane and the 
bulk, and moreover results in D-particles crossing the brane boundaries. These D-particle 
defects can even represent compactified black holes from a four-dimensional view point, with 
the extra dimensions being wrapped up appropriately in Planckian size compactifications. 
One may then encounter a situation in which D-particle point-like space time defects from the 
higher-dimensional bulk space time cross the three brane (where ordinary matter resides), 
a radically different picture from virtual excitations in a vacuum. 

In ref. js] we have discussed the details of dynamical formation of horizons on the 
brane world (in the context of (Liouville) string theory), as a result of the encounter of 
brane matter with the crossing D-particle defect. Schematically, ordinary string matter on 
the brane creates - through back reaction (recoil) effects due to scattering off D-particles- 
sufficient distortion of space-time for dynamical horizons, surrounding the defect, to appear. 
The appearance of horizons in this way looks - from the point of view of a four dimensional 
observer - as a dynamical "flashing on and off of a black hole" , coming from the "vacuum" . 
Using (weak) positive energy conditions, we have proven in [6] that such configurations with 
horizons are unstable. The life time of such objects is of the order of the Planck time, 
since this is the time uncertainty for the defect to cross the brane world and interact with 
stringy matter excitations. Once horizons form there is entropy production and through this 
irreversibility and decoherence. Consequently such stringy black hole defects are therefore 
not equivalent to ordinary virtual particles in flat space-time field theories or in attempts 
to discuss effective local quantum gravity approaches from the point of view of decoherence 
(as those mentioned in 38]). 

The presence of dynamical horizons is a real effect of the ground state of quantum gravity 
(at least in such Liouville-string approaches to QG), which implies "real" environmental 
entanglement of matter systems with (gravitational) degrees of freedom behind the horizons. 
This leads to the problem of loss of information for particles propagating outside the horizon, 
and as such can lead to microscopic time irreversibility a la Wald 39j] , and consequent CPT 
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violation and QG-induced decoherence. There is then a consequent non-unitary evolution 
of particles outside the horizon. Somewhat general arguments (even in flat space-times but 
with a boundary) have been put forward in the literature 40j to justify this point of view. 
The general message of the non-unitary evolution has then been extracted and codified with 

H | over two decades s|, 3, 33 ] to describe 
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phenomenological Lindblad master equations 
particles evolving in space-time foam. 

We cannot, of course, advocate at this stage that this (non-critical, Liouville) string 
approach, or similar, is the only consistent approach to quantum gravity. Hence we by 
no means exclude the validity of the local effective approach to QG, involving only virtual 



gravitons; in such cases there might not be any induced decoherence c 



for reasons stated 



above. It is therefore a challenging experimental issue to seek such decoherence effects 
induced by quantum gravity, which would definitely discriminate between several models of 
quantum gravity. 

Moreover, there is another interesting possibility regarding neutrinos. As pointed out 



recently in [20], the tiny mass differences between neutrino flavours may themselves (in 
part) be the result of a CPT violating quantum-gravity background. The phenomenon, if 
d be the generalisation of the celebrated Mikheyev-Smirnov-Wolfenstein (MSW) 



true, wou 
effect 41 



42] . The latter arises from effective mass differences between the various neutrino 



flavours, as a result of different type of interactions of the various flavours with matter within 
the context of the Standard Model. The phenomenon has been generalised to randomly 
fluctuating media 43] , which are of relevance to solar and nuclear reactor /3-decays neutrinos. 
This stochastic MSW effect will be more relevant for us, since we consider space-time foam, as 



a random medium which induces flavour-sensitive mass differences. If we can extrapolate |4 Z 

n 

semi-classical results on black-hole evaporation, in both general relativity [45] and string 
theory 46] to the quantum gravity foamy ground state (assuming it exists and characterizes 
the ground state of some (stochastic) quantum-gravity models, it follows that microscopic 
black holes which are near extremal (and therefore electrically charged) would evaporate 
significantly less, compared with their neutral counterparts. Thus, we may assume 20L l44j . 
that near extremal black holes in the foam would "live" longer, and as a result they would 
have more time to interact with ordinary matter, such as neutrinos. Such charged black 
holes would therefore constitute the dominant source of charge fluctuations in the foam 
that could be responsible for foam-induced neutrino mass differences according to the idea 



S 



proposed in [20]. Indeed, the emitted electrons from such black holes, which as stated above 
are emitted preferentially, compared to muons or other charged particles, as they are the 
lightest, would then have more time to interact (via coherent standard model interactions) 



with the electron-neutrino currents, as opposed to muon neutrinos 
flavour bias of the foam medium, which c ou 
gravitational analogue" of the M SW effect 



41 



d then be viewed 



20 



his would create a 



44J as the "quantum- 



42j in ordinary media (where, again, one has 



only electrons, since the muons would decay quickly). In this sense, the quantum gravity 
medium would be responsible for generating effective neutrino mass differences {^(J. Since 
the charged-black holes lead to a stochastically fluctuating medium, we will consider the 
formalism of the MSW effect for stochastically fluctuating media 43(, where the density of 
electrons would be replaced by the density of charged black hole/anti black hole pairs. 



A. Quantum-Gravitational MSW effect and induced decoherence 

For simplicity, we will give theoretical details for the case of two generations of neutrinos 
and v T with mass eigenvalues m,\ and m 2 . We take the effective Hamiltonian to be of the 
form 

H cS = H + nl h (r)Hj, (2) 

where Hi is a 2 x 2 matrix whose entries depend on the interaction of the foam and neutrinos 
and H is the free Hamiltonian. For the purposes of this paper we take this matrix to be 
diagonal in flavour space. Although we leave the entries as general constants, a Ui , we expect 
them to be of the form a GNTi^ h (r); so we write Hj as 

°y <3) 

\ a VT j 

where the foam medium is assumed to be described by Gaussian random variables [20] . We 
take the average number of foam particles, (n^ h (t)) = n (a constant), and (n^(t)n^(t')) ~ 
VL 2 n^5{t — t'). Following [43j we can deduce the modified time evolution of the density matrix 



as 

jL(p> = _ i[H + noHh {p)] _ Q%l[ Hl , [Hj, (p))) (4) 

where (...) represents the average over the random variables of the foam. The double com- 
mutator is the CPT violating term since although it is CP symmetric it induces time- 



irreversibility. It is also important to note that A2 here is of the Markovian-Liouville- 

Lindblad form for a self-adjoint operator. This is as an appropriate form for decoherence for 

environments about which we have little a priori knowledge. In the CPT violating term we 

can require the density fluctuation parameter to be different for the anti-particle sector from 

that for the particle sector, i.e. Cl ^ fl, while keeping (nl h (t)) = no the same in both sectors. 

Physically this means that neutrinos and antineutrinos with the same momenta, and hence 

interacting with the same amount of foam particles on average, will evolve differently; this 

is a result of CPT violation. 

We can expand the Hamiltonian and the density operator in terms of the Pauli spin 

matrices s M (with 4f = 1 2 the 2x2 identity matrix) as follows 

3 3 
H eS = ^2(h fi + n ti tl )-£, p = ^2 Pu ^-. (5) 

/x=0 u=0 

(where H ef[ = H + n Hj). We find that 

K = — ijfc— ^ + — a*— ^ ( 6 ) 



and 



where A; is the neutrino energy scale. 

Adopting the Lindblad decoherence approach [h]] described in 44] one can arrive to 
expressions for the pure states representing as given by 

<p)^ = il 2 + sin(2fl)| + cos(2£)| (8) 
and the corresponding state for v T as 

(p)^ = il 2 -sin(20)|-cos(2^)|. (9) 
If (p) (0) = (p) then the probability P v ^ Vt (t) of the transition — > v T is given by 

P^At)=Tr((p)(t)(p} M ). (10) 

In order to study decoherence we will calculate the eigenvectors e W and corresponding 

eigenvalues A» of £ to leading order in fl 2 . In terms of auxiliary variables U and W where 

2 2 

yyi 777. 

W = (a„„ - o„ T ) cos (20) + 1 2fc 2 (11) 
10 



and 

W= ( 

it is straightforward to show that 

~e^(i) ~ 



a, v - a„ T ) sin (20) 



W 



r {2) 



u_ . yu 2 + w 2 



... . w V^ 2 + w 2 

e (3) ~ , i , 1 

w w 



and 



Ai ~ -ft 2 ( W cos (20) - U sin (20)) ; 



Q 2 

A 2 ~ -WW + W 2 - — (W 2 + W 2 + (W cos (20) + W sin (20)) 2 ) 

Q2 

A 3 ~ iVU 2 + W 2 - — (W 2 + W 2 + (W cos (20) + W sin (20)) 2 ) . 
The vector ~~p (0) can be decomposed |44j as 

(0) = bi^w + b 2 ~? (2) + b 2 ~e > < s > 



with 



and 



Hence 



bi 



ZY 2 cos (20) +^Wsin (20) 
= U 2 + W 2 

W 2 cos (20) - UW sin (20) 



2(W 2 + W 2 ) 

p(t) = ~ (b 1 e Xlt ~t^ .~t + b 2 ~? {2) .1* + b 2 ~e*( 3 > .~t + 1 2 ) . 



(12) 



(13) 



(14) 



(15) 

(16) 
(17) 
(18) 



From this one can obtain from a standard analysis 22j, |23|, |25|, |43j the following expression 
for the neutrino transition probability <-> v T in this case, to leading order in the small 
parameter Q 2 <C 1: 



P, 



I + e -Aa^ t (l+#(eos(40)-i)) sin(t Vf) sin 2 (20)Aa 2 T ft 2 A 



2 , 3sin 2 (26l)A 2 2 

12 



4r 5 / 2 



p3/2 



e -Aa^n^(l+^(co S (40)-l)) cos(f Vf) sin 2 (2 ^) 



A 



12 



2r 



A a2 T n2 t Af 2 sin2(29) ( Aa ^ + cOs(26>)Al 2 ; 



2r 



(19) 
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a, 



where T = (Aa^ T cos(2fl) + A 12 ) 2 + Aa 2 T sin 2 (20) , A 12 = ^ and Aa MT _ ,,,,„ 

From (fl9j) we easily conclude that the exponents of the damping factors due to the 
stochastic-medium-induced decoherence, are of the generic form, for t = L, with L the 
oscillation length (in units of c = 1): 

exponent ~ -Aa 2 ^ 2 */ (0) ; f(9) = 1 + ^(cos(40) - 1) , or ^^i 26 ) 

that is proportional to the stochastic fluctuations of the density of the medium. The reader 
should note at this stage that, in the limit A12 — > 0, which could characterize the situation 
in [2^], where the space-time foam effects on the induced neutrino mass difference are the 
dominant ones, the damping factor is of the form exponent grav j tationa i msw ~ 

-fi 2 (Aa MT ) 2 L , 

with the precise value of the mixing angle 9 not affecting the leading order of the various 
exponents. However, in that follows from ffl9l) . the overall oscillation probabil- 

ity is suppressed by factors proportional to A 2 2 , and, hence, the stochastic gravitational 
MSW effect [20|], although in principle capable of inducing mass differences for neutrinos, 
however does not suffice to produce the bulk of the oscillation probability, which is thus 
attributed to conventional flavour physics. The damping exponent should then be indepen- 
dent of the mixing angle for consistency. Indeed, we find the purely gravitational MSW to 
give exponent gravitationalMSW oc Q 2 A 2 L which is independent of 9. However, this stochastic 
gravitational MSW effect, although capable of inducing neutrino mass differences, gives an 
oscillation probability which is suppressed by factors proportional to A 2 2 . Hence, the bulk 
of the oscillation is due to conventional flavour physics. 

After this theoretical discussion we now proceed to give a brief description of the most 
important phenomenological consequences of such a scenario involving decoherence. These 
can help in imposing stringent constraints on the percentage of the neutrino mass difference 
that could be due to the quantum-gravity medium. For simplicity we restrict ourselves to 
two generations, which suffices for a demonstration of the important generic properties of 
decoherence. The extension to three generations is straightforward, albeit mathematically 



more complex 25]. 



We note here that, for gravitationally-induced MSW effects (due to, say, black-hole foam 



models as in 



20 



Aa^ oc G N n 

with Gn = 1/Mp, M P ~ 10 19 GeV, the four-dimensional Planck scale, and in the case of the 

12 



gravitational MSW-like effect 20] no represents the density of charge black hole/anti-black 
hole pairs. This gravitational coupling replaces the weak interaction Fermi coupling constant 
Gp in the conventional MSW effect. This is the case that is relevant for this work. In such 
a situation the density fluctuations Q 2 can be assumed small compared to other quantities 
present in the above formulae, and an expansion to leading order in Q 2 is appropriate. 



B. Stochastic fluctuations of Space-Time metric backgrounds 

There are other models of stochastic space-time foam also inducing decoherence, for 



instance the ones discussed in 



25 



44j . in which one averages over random (Gaussian) fluc- 
tuations of the background space-time metric over which the neutrino propagates. In such 
an approach, one considers merely the Hamiltonian of the neutrino in a stochastic metric 
background. The stochastic fluctuations of the metric would then pertain to the Hamiltonian 
(commutator) part of the density-matrix evolution. In parallel, of course, one should also 
consider environmental decoherence-interactions of Lindblad (or other) type, which would 
co-exist with the decoherence effects due to the stochastic metric fluctuations in the Hamil- 
tonian. For definiteness in what follows we restrict ourselves only to the Hamiltonian part, 
with the aim of demonstrating clearly the pertinent effect and study their difference from 
Lindblad decoherence. 

In this case, one obtains transition probabilities with exponential damping factors in 
front of the oscillatory terms, but now the scaling with the oscillation length (time) is 



quadratic 



44|, consistent with time reversal invariance of the neutrino Hamiltonian. For 
instance, for the two generation case, which suffices for our qualitative purposes in this work, 
we may consider stochastically fluctuating space-times with metrics fluctuating along the 



direction of motion (for simplicity) 44] 



ar 



-(ai + l) 2 + al -a 3 (ai + 1) + a 2 (a 4 + 1) 
-a 3 (ai + 1) + a 2 (a 4 + 1) -a| + (a 4 + l) 2 



with random variables (a*) = and (a* a?) = S^ai. 



Two generation Dirac neutrinos, then, which are considered for definiteness in 44j (one 
would obtain similar results, as far as decoherence effects are concerned in the Majorana 
case), with an MSW interaction V (of unspecified origin, which thus could be a space-time 
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foam effect) yield the following oscillation probability from an initial state of flavour 1 to 2 
is 

Prob(l -> 2) = UvUqUuUue'te-"** (20) 
hi 

where the time dependent part is 

with U the mixing matrix, which, in the two-flavour-dominance scenario we are working on 
here for the sake of brevity, can be parametrised by a mixing angle 9: 

cos 9 sin 9 
— sin 9 cos 9 

Since the {a,} are assumed to be independent Gaussian variables, the pertinent covariance 
matrix H has the diagonal form 

/ J- \ 



u 



(21) 



0^00 



(22) 



0^0 

V o o o 

with o"j > 0. The calculation of transition probabilities requires the evaluation of the follow- 
ing average over the stochastic space-time fluctuations dj 

,det H 



The result is 



44 1 



) = d a exp(— a • S • a)e 



7T- 



(23) 
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X 
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2p 



,„ , 2 V cos 2£f m? — m~ ) ,^ „ „ ^ \ » 

(9ffi +<X 2 +<T 3 +CT4 )+ H 21 (12<xi +2cr 2 -2<r 3 ) ) * 2 



(24) 



where again is the neutrino energy scale, <Ji ,i = 1, . . .4 parametrise appropriately the 

yfc ITI < 1, and A; 2 > 



stochastic fluctuations of the metric in the model of [44], T 

m 2 , mf, and 



15 7 ^ 



4=2 (m? + T(l + cos2#)(m 2 - m|) + T 2 (m 2 - m|) sin 2 20) 



1 



(m 2 + T(l - cos2fl)(m 2 - m\) - T 2 (m 2 - m\) sin 2 20) 



(25) 
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Note that the metric fluctuations-dj induced modifications of the oscillation period, as well as 



exponential e - '- "-'* time-reversal invariant damping factors 44], in contrast to the Lindblad 
decoherence, in which the damping was of the form e - '"")*. At first thought, one may 
attribute this feature to the fact that, in this approach, only the Hamiltonian terms are 
taken into account (in a stochastically fluctuating metric background), and as such time 
reversal invariance t — > — t is not broken explicitly. But there is of course decoherence, and 
the associated damping. 



However, upon closer inspection things are not as simple. As shown in [36], the power 
of the time variable t in the associated damping is crucially dependent on the type of the 
distribution characterising the gravitational fluctuations. In terms of our D-particle-recoil 
induced stochastic model of space-time foam 0, Q, 44], such distributions characterise 
the ensemble of velocities of the gas of D-particles involved in the foam, with which the 
neutrinos interact, and which in turn affects the induced metric fluctuations, as explained 
above. Assume, for instance, that the distribution of the (induced) metric fluctuations are 
of Cauchy-Lorentz type, which could be induced by a distribution of D-particle velocities 



in a D-foam model of the type considered in 



3- 



Such a distribution has undefined mean 



and variance as well as undefined or infinite higher moments. Assuming for concreteness a 
case with zero mean ( in the sense of principal values), the pertinent distribution function 



is taken to be 



36|: 



m = ^ (26) 

with £ > the characteristic scale parameter of the Cauchy-Lorentz distribution. In that 
case, the pertinent statistical average of the associated oscillation probability for a two- 
flavour neutrino problem would be given, to leading order in an expansion in powers of 
rrii/k, of interest to us here, by: 

(e *K- 2 )t )cL „ exp {tktA _ £fct|A|) , A = m '~ 2 m ' ( 27 ) 

From the above equation we do observe a linear damping, similar to the Lindblad environ- 
ment case. The important point to notice is that in this case, the damping exponent is of 
order 

exp (-nl L t) , with fi* L ~ g 1 ™*"™ 11 (28) 



The reader should compare the linear power of the (small) quantity \m\ — m\\/k entering the 
damping exponent (1281) in the Cauchy-Lorentz stochastic foam model to the quadratic power 
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of that quantity entering the Gaussian model of foam (1241) . where the pertinent exponents 
are proportional to factors of the form (m 2 — vn^) 2 jk 2 and hence much smaller, provided 
of course the parameters a 2 and £ are of similar order. At this stage, these are treated as 
phenomenological parameters, since their order depends on the details of the underlying 
model. For instance, in the case of the D-particle foam model 0, , such an information 
depends on the dynamics of the gas of bulk D-particles, which probably is an issue that can 
only be resolved within a microscopic M-theory model. 



At this stage we would like to draw the reader's attention to a possible interpretation [29] 
of the Lindblad-type exponential damping ( 128]) of the Cauchy-Lorentz distribution with 1/E 
dependence as a neutrino decay, 

exp(-fi^ L t) = exp(-t/r lab ) = exp(-tm u J Er TCSt ), (29) 

from which we can get a lower limit on the (unstable) neutrinos life times. We shall come 
back to this issue in the concluding section of the article, when we provide the relevant 
experimental bounds. 

It must be stressed at this point, before closing this subsection, that the above consid- 
erations, especially the ones concerning the form and order of the decoherence damping 
factors, of interest to us in this work, although derived in a two-dimensional toy-model 
(considering metric deformations primarily along the direction of motion of the neutrino 
probe), nevertheless are valid qualitatively in a full fledged four space time dimensional 
model. This has been demonstrated in [36J], where realistic models of neutrinos propagating 
in four-dimensional, stochastically fluctuating, space-time backgrounds have been considered 
in detail, with results similar to the ones considered in 12J and reviewed above. 



C. Mimicking Decoherence via Conventional Uncertainties in Neutrino Energy 
and Oscillation Length 

A few remarks are now in order regarding the similarity of this latter type of decoherence 
(1241) with the one mimicked [^tJ by ordinary uncertainties in neutrino experiments over the 
precise energy E of the beam (and in some cases over the oscillation length L). Indeed, 
consider the Gaussian average of a generic neutrino oscillation probability over the L/E 
dependence (P) = dxP(x) ^^^ e~ ^ 2 , with I = (x) and a = \J JJx — (x)) 2 ), x = j^, 
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and assuming the independence of L and E, which allows to write (L/E) = (L)/{E). A 



pessimistic and an optimistic upper bound for a are given by 



47) 



pessimistic: a ~ Ax = A^ < AL|f§ \ L={L) ^ E={E) + &E\§e\l=(L),e=(e) 



(L) ( AL , AE \ 

me) y (l) <£> ; 



. optimistic: v<i%\/(^y+(^ 2 
Then, it is easy to arrive at the expression 471 ] 

(Pa 13) = $a/3 



(E) 



n 11 

2 E E Re (KaUp a U ab U; b ) (l - co S (2lAml b )e- 2 ^ Am ^ 2 

a=l b=l,a<b 
n n 



E Im toaUpaUcbUfr) sin(2Mm^)e- 2 ^ A ^ 



, ) 

ll " \^ aa^ pa^ ao^ pbj u'"v i "-^"" , a6/'' 

a=l b=l a<fe 



with (30) 

with [/ the appropriate mixing matrix. Notice the a 2 damping factor of neutrino oscillation 
probabilities, which has the similar form in terms of the oscillation length dependence (L 2 
dependence) as the corresponding damping factors due to the stochasticity of the space-time 
background in ff24)) . It is noted, however, that here / has to do with the sensitivity of the 
experiment, and thus the physics is entirely different. 

In the case of space-time stochastic backgrounds, one could still have induced uncertain- 
ties in E and L, which however are of fundamental origin, and are expected to be more 
suppressed than the uncertainties due to ordinary physics, described above. Apart from 
their magnitude, their main difference from the uncertainties in (I30p has to do with the 
specific dependence of the corresponding a 2 in that case on both E and L. For generic 
space-time foam models it is expected that an uncertainty in E or L due to the "fuzziness" 
of space time at a fundamental (Planckian) level will increase with the energy of the probe, 
SE/E,SL/L oc (E/Mp) a , a > 0, since the higher the energy the bigger the disturbance 
(and hence back reaction)on the space time medium. In contrast, ordinary matter effects 



decrease with the energy of the probe [47J, |48|. 
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III. PREVIOUS DECOHERENT FITS WITH EXISTING NEUTRINO DATA 
AND PHYSICAL INTERPRETATION 



The first complete phenomenological attempt to fit decoherence models to atmospheric 
neutrino data was done in [ijj], where for simplicity a two-generation neutrino model with 
completely positive Lindblad decoherence, characterised by a single parameter 7, and lead- 
ing to exponential damping with time of the relevant oscillatory terms in the respective 
oscillation probabilities, was considered. 

Various dependencies on the energy E of the neutrino probes have been assumed, in 
a phenomenological fashion, for the Lindblad decoherence coefficient 7 = 7L n b (gIv)" > 
with n = 0,2,-1. The sensitivities in the work of [19] from atmospheric neutrinos (plus 



261 ]) at 90% C.L. can be summarised by the following bounds on the 



accelerator data 
parameter 7Lnt> : 

7Lnb < 0.4 x 10~ 22 GeV , n = 

7Lnb < 0.9 x 10~ 27 GeV , n = 2 

7Lnb < 0.7 x 10~ 21 GeV , n = -1 (31) 

Recently {27], updated values on these parameters, refered to 95% C.L., have been provided 
by means of combining solar-neutrino and KamLand data 

7Lnb < 0.67 x 10~ 24 GeV , n = 
7 Lnb < 0.47 x 10~ 20 GeV , n = 2 

7Lnb < 0.78 x 10~ 26 GeV , n = -1 (32) 

It should be remarked that all these bounds should be taken with a grain of salt, since there 
is no guarantee that in a theory of quantum gravity 7L n b should be the same in all channels, 
or that the functional dependence of the decoherence coefficients 7 on the probe's energy E 
follows a simple power law. Complicated functional dependencies ^y(E) might be present. 

We shall come back to these bounds in the discussion section of the article, when we 
compare the potential of upcoming neutrino data, with energies of order of tens of GeV, 
from CNGS facility. We also investigate the sensitivity of the experiments at J-PARC. The 
J-PARC beam operates at rather lower energies comparing to CNGS, however the fact that 
the maximum of oscillation in the spectrum will be measured by T2K experiment, allows 
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FIG. 1: Left : Ratio of the observed v e spectrum to the expectation versus Lq/E for our decoherence 
model. The dots correspond to KamLand data. Right: Decoherence fit. The dots correspond to 
SK data. 



to achieve a remarkable sensitivity to those dapmping exponents with low power energy 
dependence, as compared with [19j, |27|. The CNGS is very sensitive to the E 2 dependent 
case despite the fact that the spectrum of atmospheric neutrinos used in JjJ spans a wide 
range of energy which extends to 100-1000 GeV. The E 2 dependence, for instance, could 
characterise Cauchy-Lorentz stochastic models of space time foam. 

In 25j a three generation Lindblad decoherence model of neutrinos has been compared 



against all available at the time experimental data, taking into account the recent results 



from KamLand experiment 



341 ] indicating spectral distortions. 



The results are summarized in Fig. [TJ which demonstrates the agreement (left) of the 
model with the KamLand spectral distortion data and the best fit (right) for the 
Lindblad decoherence model used in ref. 

The best fit has the feature that only some of the oscillation terms in the three generation 
probability formula have non trivial damping factors, with their exponents being independent 
of the oscillation length, specifically 25|. If we denote those non trivial exponents as T> ■ L, 



we obtain from the best fit of 



25|: 



V 



1.3 ■ 10~ 2 



(33) 
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in units of 1/km with L = t the oscillation length. The 1/L-behaviour of T>n, implies, as 
we mentioned, oscillation-length independent Lindblad exponents. 

In 25] an analysis of the two types of the theoretical models of space-time foam, discussed 



in section 2, has been performed in the light of the result of the fit (1531) . The conclusion was 
that the model of the stochastically fluctuating media (1191) (extended appropriately to three 
generations [25|, so as to be used for comparison with the real data) cannot provide the full 
explanation for the fit, for the following reason: if the decoherent result of the fit (1331) was 
exclusively due to this model, then the pertinent decoherent coefficient in that case, for, say, 
the KamLand experiment with an L ~ 180 Km, would be \V\ = VL 2 G 2 N n 2 ~ 2.84- 10~ 21 GeV 
(note that the mixing angle part does not affect the order of the exponent). Smaller values 
are found for longer L, such as in atmospheric neutrino experiments. The independence of 
the relevant damping exponent from the oscillation length, then, as required by fl33l) may be 
understood as follows in this context: In the spirit of 20], the quantity Gj^uq = £^r-, where 
£ <C 1 parametrises the contributions of the foam to the induced neutrino mass differences, 
according to our discussion above. Hence, the damping exponent becomes in this case 
( 2 fi 2 (Am 2 ) 2 • L/E 2 . Thus, for oscillation lengths L we have L~ x ~ Am 2 /E, and one is left 
with the following estimate for the dimensionless quantity £ 2 Am 2 Q 2 /E ~ 1.3 ■ 10~ 2 . This 
implies that the quantity Q 2 is proportional to the probe energy E. In principle, this is not 
an unreasonable result, and it is in the spirit of 20(], since back reaction effects onto space 
time, which affect the stochastic fluctuations Q 2 , are expected to increase with the probe 
energy E. However, due to the smallness of the quantity Am 2 /E, for energies of the order 
of GeV, and Am 2 ~ 10 -3 eV 2 , we conclude (taking into account that £ <C 1) that Q 2 in this 
case would be unrealistically large for a quantum-gravity effect in the model. 

We remark at this point that, in such a model, we can in principle bound independently 
the Q and no parameters by looking at the modifications induced by the medium in the 
arguments of the oscillatory functions of the probability ( fl9l) . that is the period of oscillation. 
Unfortunately this is too small to be detected in the above example, for which Aa eM <C A12. 

The second model ( 1241 of stochastic space time can also be confronted with the data, 
since in that case (1331) would imply for the pertinent damping exponent 
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K - mn) 2 ; , 2Kcos26>(m 2 — ml) . A , 

1 J -(9(7! +(72 + ^3 + <7 4 ) + H -(12(71 + 2(X 2 - 2<T 3 ) * 



~ 1.3 • 10" 2 . (34) 

Ignoring subleading MSW effects V, for simplicity, and considering oscillation lengths 
t = L ~ (mt-rn 1 ) ' we then observe that the independence of the length L result of the 
experimental fit, found above, may be interpreted, in this bounding the stochastic 

fluctuations of the metric to %a\ + o~2 + 03 + 04 ~ 1.3. ■ 10~ 2 . Again, this is too large to 
be a quantum gravity effect, which means that the L 2 contributions to the damping due 



to stochastic fluctuations of the metric, as in the model of 
exclusive explanation of the fit. 



441 ] above (1241) . cannot be the 



The analysis of 25| also demonstrated that, at least as far as an order of magnitude of 
the effect is concerned, a reasonable explanation of the order of the damping exponent (1331) , 
is provided by Gaussian-type energy fluctuations, due to ordinary physics effects, leading 
to decoherence-like damping of oscillation probabilities of the form (1301) . The order of these 
fluctuations, consistent with the independence of the damping exponent on L (irrespective 



of the power of L), is 

AE 



E 



1.6 - (35) 



if one assumes that this is the principal reason for the result of the fit. 

However, not even this can be the end of the story, given that the result (|33|) pertains 
only to some of the oscillation terms and not all of them, which would be the case expected 
for the ordinary physics uncertainties ( 1301) . The fact that the best fit model includes terms 
which are not suppressed at all calls for a more radical explanation of the fit result, and the 
issue is still wide open. 

It is interesting, however, that the current neutrino data can already impose stringent 
constraints on quantum gravity models, and exclude some of them from being the exclusive 
source of decoherence, as we have discussed above. Of course, this is not a definite conclusion 
because one cannot exclude the possibility of other classes of theoretical models of quantum 
gravity, which could escape these constraints. At present, however, we are not aware of any 
such theory. We would like now to revisit the above constraints in upcoming neutrino data 
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from the experiments at the CNGS and J-PARC facilities. This is discussed in the next 
sections. 



IV. THE COMBINED FIT TO QUANTUM-GRAVITY DECOHERENCE SIGNA- 
TURES 

In the previous sections we have discussed several theoretical models of quantum-gravity- 
induced decoherence independently, assuming each time only one dominant type of deco- 
herence: (i) Lindblad-type, through the representation of the quantum gravity space-time 
foam as a stochastic medium, (|T9|) . (ii) stochastically fluctuating space-time backgrounds, 
fl24"|) . and (iii) induced decoherence-like evolution, as a result of uncertainties in the energy 
and/or oscillation lengths of the neutrinos, (|30|) . 

The various types of decoherence can be mainly distinguished by the form of their ex- 
ponential damping factor, as far as the power of the oscillation length L in the exponent is 



concerned, and the associated energy dependence [49j. Model independent data fits should 
combine, in general, the various types of decoherence-deformed oscillations, given that dom- 
inance of one or the other type may not be necessarily a feature of a quantum-gravity model. 

It is the purpose of this section, and one of the main objectives of this work, to establish 
the limit of sensitivity of CNGS and J-PARC beams, in a model-independent way, to a 
simple parametrization of the above effects, combined in a single model for oscillations 
between flavours a, b — l...n of the form: 

{Pap) — &ap — 
n n 

2 E E Re KaUp a U ab U; b ) (l - cos(2£Am^)e-^ i2 

a=l b=l,a<b 
n n 

E Im {U: a Up a U ah U; b ) sm(2eAml b )e-^ L2 

a=l b=l,a<b 

with i = A (36) 
where L ~ t (in units of c = 1) is the oscillation length. In general one may parametrize 



the damping exponents by polynomials in L 



49] of any degree, but parametrisations of 



degrees higher than 2 are not favoured by the class of quantum-gravity decoherence models 



considered in the literature so far [13, 



441 ]. and reviewed above. 
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From (T19]) . (I24I) . we observe that (136!) is oversimplified in that it ignores possible modi- 
fications of the oscillation period, which do exist in various microscopic models as a result 
of the decoherence or stochastic-medium effects. A complete theoretical treatment requires 
solving the evolution equations for the reduced density matrix of neutrinos in a combined 
situation involving simultaneously stochastic fluctuations of the background space-time met- 
ric and interactions (of Lindblad type) with a stochastically-fluctuating quantum-space-time 
medium. This will be left for future work. However, for our purposes in the current article, 
we note that it is a reasonable assumption that such modifications to the oscillation period 
are suppressed as compared with the ordinary oscillation terms, and as such the dominant, 
model-independent, terms appear to be only the exponents of the damping factors. Concern- 
ing the latter, we also observe from fTiTJl) that, in general, there are slight differences among 
the various exponents accompanying the oscillation terms in stochastic-medium models, 
which however are all of the same order of magnitude, and hence the error one makes in 
assuming the simplifying two-parameter (gi, ^2) damping decoherence form ( |36l) is negligible. 

For our phenomenological purposes in this work, therefore, the only important point to 
notice is that the parameters qi,i = 1,2 may be in general energy dependent, expressing 
back-reaction effects of the (neutrino) matter onto the fluctuating space-time. Following 

n □ 

earlier treatments and theoretical quantum-gravity-decoherence models [13J, [19( we shall 
consider the following three cases of generic energy dependence of the decoherence coefficients 
qui = 1,2: 

ft, i = 1,2 oc E n , n = -1,0,2 (37) 

where the reader should have in mind that in each case the pertinent decoherence coefficient 
has the appropriate units, as being a dimensionful quantity. 

For our studies we use two sets of the one and two parametric models covering the main 
variety of phenomenologies for quantum gravity induced decoherence phenomena described 
by the expression (|36l) . The first set of the models under consideration concerns the pres- 
ence of linear Lindblad-type mapping operator in the equation for the evolution of the 

density matrix for the pure neutrino quantum states ,Q,y, y,y,Q- 

The oscillation 

probabilities corrected for the decoherence effects with different energy dependence in the 
exponentials read 
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no neutrino-energy dependence 
= l sin 2 (29 23 ) 



, 2.54Am 2 
1 - exp(-5 ■ 10 7 L) cos ( — L 



(38) 



inversely proportional to the neutrino energy (e.g. the case of Cauchy-Lorentz type of 



stochastic foam 



36i, (BED 
1 



sin 2 (20 23 ) 



-2.5472 L /2.54Am 2 
1 — exp( ) cos I L 



E 



E 



(39) 



proportional to the neutrino energy squared 
1 



Pu 



sin 2 (2# 23 ) 



1 - exp(-5 ■ 10 27 l2 E 2 L) cos ( 2 ' 54 ^ m ' L 

E 



(40) 



where 70, 7 2 _ 1 and 72 are measured in eV, eV 2 length and eV -1 respectively the mass 
square difference Am 2 , is measured in eV 2 , the energy E, is measured in GeV; and the path, 
L, is measured in km. 

The second set of the models concerns the gravitational MSW stochastic effect ffl9l) with 
linear and quadratic (I24"I) time dependent fluctuations of space-time foam described by 



Pu 



1 



exp(— K\ 



cos 2 (26 : 



23 J 



— - exp(— K2) cos 



r 2 r-v 2 

where the exponential damping factors are chosen as 
• no energy dependence, with linear 



2. 54 Am 2 
E 



L sin 2 (20 23 ), (41) 



ki = 5 • 10VLsin 2 (2fl); k 2 = 5 • 10 9 « 2 L(1 + 0.25(cos(4#) - 1)) (42) 



9^2; 



quadratic 



«i = 2.5 ■ 10 i9 a 2 L 2 sin 2 (2#); k 2 = 2.5 ■ 10 19 a 2 L 2 (l + 0.25(cos(4#) - 1)) (43) 



and combined time evolution 



Kl = (5 • 10 9 7l 2 L + 2.5 ■ 10 19 7 2 L 2 ) sin 2 (29); 
K2 = (5 - 10 9 7 2 L + 2.5 ■ 10 19 7 2 L 2 )(1 + 0.25(cos(4#) - 1)) 



(44) 
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• proportional to the neutrino energy, with linear 

Kl = 5 • 10 18 f3 2 EL sin 2 (26); k 2 = 5 ■ 10 18 (3 2 EL(1 + 0.25(cos(4fl) - 1)) (45) 

quadratic 

Ki = 2.5 • 10 28 [3 2 EL 2 sin 2 (2#); k 2 = 2.5 • W 28 [3 2 EL 2 (1 + 0.25(cos(4#) - 1)) (46) 

and combined time evolution 

Kl = (5 • 10 18 7 f EL + 2.5 • 10 28 72 2 £L 2 ) sin 2 (20); 
kjj = (5 • 10 18 7f EL + 2.5 • 10 28 72 2 £L 2 )(1 + 0.25(cos(4#) - 1)) (47) 

• proportional to the neutrino energy squared, with linear time evolution 

Ki = 5 • 10 27 (3 2 E 2 L sin 2 (2#); k 2 = 5 • 10 27 (3 2 E 2 L(1 + 0.25(cos(4#) - 1)) (48) 

The energy and the path length in (I42p - (j4"8"|) are measured in GeV and km respectively, 
while the parameters in damping exponentials are given in eV in respective power (see 
Tabled] for details). 

V. SENSITIVITY OF CNGS AND J-PARC BEAMS TO QUANTUM-GRAVITY 
DECOHERENCE 

In this section we study the expected sensitivity of the CNGS and J-PARC beams to the 
quantum gravitational decoherence phenomena described by ( |38l) -( l4*8l) . considering them as 
subdominant contributions to the atmospheric oscillations effects. 

Both CNGS and J-PARC are conventional neutrino beams where neutrinos are produced 
by the decay of secondary particles (pions and kaons) obtained from the collision of the 
primary proton on a graphite target. For the CNGS beam, the protons come from the 
CERN-SPS facility with a momentum of 400 GeV/c whereas in the case of the J-PARC [50| 
the protons are produced in Tokay (Japan) and have a momentum of 40 GeV/c. The 
expected number of protons on target per year at the nominal intensity is 4.5 x 10 19 and 
1 x 10 21 respectively for the CNGS and J-PARC beam and the envisaged run length is 5 
years in both cases. 
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Both beams will be used for long baseline neutrino experiments which, starting from a 
beam, will search for neutrino oscillations. The OPERA experiment will measure neutrino 
events on the CNGS beam using a 2 kton detector which relies on the photographic emulsion 
technique, located at a baseline of 732 km; the first neutrino events were observed in August 
2006 [5lt 

The T2K experiment will use the J-PARC beam measuring neutrino events with the 



Super-Kamiokande 52| detector (a water cerenkov detector with an active volume of 22.5 
kton) at a baseline of 295 km. 

Although CNGS beam designed in a way to be optimised for the — *> v T oscillation 
searches through the detection of r lepton production in a pure beam there is also a 
possibility to measure spectrum by reconstructing /i from the charged current (CC) 
events caused by v^. Moreover, for this experiment, we can take advantage of high mean 
value for the energy of z/ M s which makes the exponential damping factors more pronounced 
for some cases described in the previous section. 

The number of /i is given by the convolution of the flux d(j) Vll /dE with the CC cross 
section on lead cr^ c (P), weigted by the — > surviving probability P u ^ Ufi , times the 
efficiency e w of muon reconstruction of a given detector: 

dN ^ = A ^P a cc (E)e (49) 

where is a normalisation factor which takes into account the target mass and the nor- 
malisation of the Vn in physical units. In our study we assumed an overall efficiency e w of 
93.5% for the OPERA experiment and of 90% for the T2K one as stated in the experiment 
proposals. 

To estimate quantitatively the sensitivity of CNGS on P v ^ Vt described by ( l38l) -( H8l) . we 
simulated the theoretical spectra of the reconstructed events for various values of damping 
parameters. Since there is no near detector at the neutrino source the overall normalisation 
of the un-oscillated neutrino flux cannot be controlled with the precision better than 20%, 
therefore such a normalisation has been taken into account in our x 2 analysis to estimate 
the expected limits of sensitivity on the damping parameters: 



x 2 



J2[x t -aP t } 2 /a 2 + (l-a) 2 /a 2 , (50) 



where, x« is the expected number of CC events contained in the i-th energy bin considering 
standard three flavour oscillation, Pj is the number of events in the i-th bin theoretically 
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expected when some decoherence parameters are considered and Ui represents the error on 
the number of events in the i-th bin. The parameter a represents the normalisation factor 
and the additional contribution (1 — a) 2 /a 2 is related to the systematic uncertainty of the 
overall neutrino flux at the source (a = 0.2). This systematic unceartanty 60j plays an 
important role in the correct estimation of the sensitivity of the experiment especially when 
the shape of CC events spectrum is not changed by the decoherence effects (i.e. exponents 
independent on energy). 



For the atmospheric parameters we used Am 2 = 2.5 • lCT 3 eV 2 and #23 = 45° 35|. The 
3 cr sensitivity on the damping parameters is found by applying a cut on the value of the \ 2 
of 9 and 11.83 respectively for 1 d.o.f and 2 d.o.f. 

As the CNGS beam is designed to observe u T , neutrinos will have a high energy with a 
mean value of about 17 GeV. This represents an advantage since it makes the exponential 
damping factors more pronounced for some cases described in the previous section. For the 
OPERA experiment the systematic uncertainty in the muon detection efficiency is negligible 
compared to the statistical uncertainties, therefore the error o"j used in our analysis (Eq. 1501) 
represents the statistical error only. 

To generate the expected neutrino spectra of the CNGS beam measured by the OPERA 
experiment we used a fast simulation algorithm described in 53j (see also Appendix A 
for details). We present in Figj2] a typical simulated spectrum of the expected number 
of \x events including the effects of decoherence (for the case of an inversely proportional 
dependence on neutrino energy) as a sub dominant suppression of the probability inferred 
from the atmospheric neutrino experiment 351 ] . 

Our results for the sensitivity of CNGS to one parametric decoherence damping expo- 
nentials in P v -+ Vt are summarised in second column of Table [B Also, for two parametric 
fits (jHJ) and flSj) , the 3 a CL sensitivity contours are presented in Fig. @] and Fig. [51 

Contrary to the OPERA experiment, the T2K experiment was designed to observe v e 
and the mean energy is much lower: the maximum of oscillation at the given baseline of 295 
km corresponds to a neutrino energy of about 600 MeV and a narrow spectra at the selected 
energy will be obtained using the so called off-axis technique [54]. The spectrum covers the 
region of the first maximum of oscillation and this is a region where the QG effects could be 
easily observed due to the small number of CC events expected in case of no QG damping 
exponents, as it can be seen in Fig [3J 
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FIG. 2: The number of reconstructed v„ CC events in OPERA as a function of the neutrino energy 
with (blue line) and without (red line with error bars) QG decoherence effect included in case of 
inversely proportional dependence on neutrino energy. 3a difference between the expected and QG 
disturbed spectra is shown. 

The neutrino production at J- PARC beam is simulated in GEANT environment, which 
takes into account the whole focusing system (horn and reflectors), target and decay tunnel 
at J-PARC. Protons are generated on target and through the decay of parent pions and kaons 
the probability of neutrino at a selected location is calculated and the spectra is obtained. 
We use the reconstructed neutrino energy for single-Cherenkov-ring muon quasi-elastic (QE) 
and non-QE events. Of course, here the energy resolution plays an important role: for this 
reason we introduced an energy smearing effect of 20% in our analysis. This value takes 
into account the different energy resolution for the two kind of events and the fact that QE 
events are the majority of the of muon neutrino events in the detector. 

Our results obtained using the same way of analysis quantified by (149]) and (150~1) for 
the sensitivity of T2K to one parametric decoherence damping exponentials in P v are 
summarised in third column of Table [H Also, for two parametric fits (jHJ) and (14 7p . the 
3 a CL sensitivity contours are presented in Fig. [6] and Fig. [7J Contrary to CNGS, the 
J-PARC facility is equipped with a near detector which measures the un-oscillated muon 
spectrum with 5% uncertanty in the absolute normalization of the overal flux. However, to be 
conservative, we obtain our results in 3rd and 4th column of Table H under the assumption 
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FIG. 3: The number of reconstructed v„ CC events in as a function of the neutrino energy with 
(blue line) and without (red line with error bars) QG decoherence effect included in case of inversely 
proportional dependence on neutrino energy. 3a difference between the expected and QG disturbed 
spectra is shown. 



of 20% uncertanty in the overal normalization of the spectrum. Since the main effect is 
related to the maximal oscillation point in the spectrum, the overall normalization is not as 
critical as in the case of CNGS fit. 

The T2K experiment yields a better limit on the damping parameters only in the case 
where the effect has no energy dependence or contains inversely proportional to the neutrino 
energy exponent, as expected given the low energy spectrum. In all the other cases, the 
dependence on the baseline disfavours the short baseline of T2K with respect to OPERA. 

Another possibility to observe the effect on the T2K neutrino beam is to select a longer 
baseline, namely to locate the detector at about 1000 km in Korea. Studies of beam up- 



grades and a large liquid Argon detector of 100 kton in Korea were carried out [55j in the 
framework of CP violation discovery. We considered this option, called T2KK, and studied 
the possibility to constrain damping parameters in this case. The proposed upgrade at 4 
MW of the beam was taken into account which results into 7 x 10 21 p.o.t. per year and a 
running time of 4 years was envisaged. The efficiency e w for the detector is assumed to be 
95% and an energy smearing of 15% is taken into account. 

Our results for the sensitivity of T2KK to one parametric decoherence damping exponen- 
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tials in P v _> Vt are summarised in fourth column of Table [B This configuration yields better 
results than the T2K experiment and results comparable to the OPERA experiment. 

VI. DISCUSSION AND OUTLOOK 



It is instructive to compare the sensitivity limits presented in Table [I] with those derived 
from the analysis of atmospheric neutrino data [19] obtained at Super-Kamiokande and 
K2K experiments. One can transform the limits on the Lindblad type operators presented 
in Table [I] to the notations ([3"Tj) of [19!] using the following transformations: 

7Lnb = 7o[GeV], n = 

7Lnb = 7 2 [GeV- 1 ] x (GeV 2 ), n = 2 



7Lnb = 7 'i[GeV 2 ]/(GeV), n 



(51) 



Lindblad-type mapping operators 


CNGS 


T2K 


T2KK 


70 [eV] ; ([GeV]) 


2 x 10" 13 ; (2 x lCT 22 ) 


2.4 X 10~ 14 ; (2.4 X 10" 


-23) 


1.7 x 10" 14 ; (1.7 x 10" 


-23\ 


7*! [eV 2 ] ; ([GeV 2 ]) 


9.7 x 1CT 4 ; (9.7 x 10" 22 ) 


3.1 x 1CT 5 ; (3.1 x 1CT 


23) 


6.5 x 10" 5 ; (6.5 x 10" 


23) 


72 [eV 1 ] ; ([GeV- 1 ]) 


4.3 x 10" 35 ; (4.3 x 1CT 26 ) 


1.7 x 10~ 32 ; (1.7 x 10" 


-23) 


3.5 x 10" 33 ; (3.5 x 10" 


-24) 


Gravitational MSW (stochastic) effects 


CNGS 


T2K 


T2KK 


o? 


4.3 x lO" 13 eV 


4.6 x 10 -14 eV 




3.5 x 10" 14 eV 




<A 


1.1 x 10~ 25 eV 2 


3.2 x 10" 26 eV 2 




6.7 x 10~ 27 eV 2 




p 2 


3.6 x 10~ 24 


5.6 x 10~ 23 




1.7 x 10~ 23 




Pi 


9.8 x 1(T 37 eV 


4 x 10" 35 eV 




3.1 x 10" 36 eV 




Pi 


8.8 x 1CT 35 eV -1 


3.5 x 10" 32 eV -1 




7.2 x 10" 33 eV -1 





TABLE I: Expected sensitivity limits at CNGS, T2K and T2KK to one parametric neutrino deco- 
herence for Lindblad type and gravitational MSW (stochastic metric fluctuation) like operators. 



30 




FIG. 4: The expected CNGS sensitivity contour at 3 a CL, with two decoherence parameters 
contributing to the combined time evolution of gravitational MSW (stochastic metric fluctuations) 
damping with no energy dependence. 



)e directly compared with the bounds fl3Tj) 



ljj (see for details (J3H)) at 95% C.L. on 



so that the numbers of Table [J in parentheses can 
and (132|) . In particular, the bound obtained in 
the Lindblad type operators with no energy dependence is close to the sensitivity estimated 
in our analysis in case of T2K and T2KK simulations. Although, the CNGS estimation is 
about an order of magnitude weaker, one should stress that the current limit is given at 99% 
C.L. under the assumption of the most conservative level of the uncertainty of the overall 



neutrino flux at the source. The bound on the inverse energy dependence given in [19] (!3T|) 
is close to the current CNGS estimates. T2K and T2KK demonstrate an improvement. In 
spite of the fact that the Super-Kamiokande data contains neutrino of energies up to ~TeV, 
the sensitivity one obtains at CNGS to the energy-squared dependent decoherence is close, 
within an order of magnitude, to the bound fl3Tl) imposed by atmospheric neutrinos and sur- 
passes T2K and T2KK sensitivity bounds by ~ 3 and ~ 2 orders of magnitude respectively. 
The much less uncertain systematics of CNGS compared to the atmospheric neutrino data 
will make the expected bound more robust as soon as the upcoming data from OPERA will 
be analysed. Moreover, our results are also competitive with the sensitivity to the same 



Lindbland operators estimated in [56] for ANTARES neutrino telescope, which is supposed 



to operate at neutrino energies much higher than CNGS and J-PARC experiments 61|. 
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FIG. 5: The expected CNGS sensitivity contour at 3 a CL, with two decoherence parameters 
contributing to the combined time evolution of gravitational MSW (stochastic metric fluctuations) 
damping which are proportional to the neutrino energy. 

Assuming that the decoherence phenomena affect all particles in the same way, which 
however is by no means certain, one might compare the results of our analysis with bounds 



obtained using the neutral kaon system 57| . The comparison could be done for the constant 



(no-energy dependence) Lindblad decoherence model. The main bound in 57| in such a 
case reads 70 < 4.1 x 10~ 12 eV, thus being about two orders of magnitude weaker than the 
sensitivity forecasted in the present paper. 

Finally, we compare the estimated sensitivity with the bounds obtained in [^j using 
solar+KamLAND data. In principle, as in the case of the neutral kaon system, a direct 
comparison is impossible, since the parameters investigated here for the — > v T channel 
need not be the same for the v e — > channel. However, again, if these parameters are 
assumed to be roughly of equal size, then one can see that the estimates of [27( (1321) . which 
win essentially over the CNGS, T2K and T2KK sensitivities only for the case of inverse 
energy dependent decoherence, which strongly favours low neutrino energies (e.g. the case 



of Cauchy-Lorentz stochastic space-time foam models of 36( (1281) . for which the current 
limit would bound, on account of (|39|) . t 
£ < 5 x 1CT 3 for neutrino-mass differences [27j |m 2 — m 2 \ = (7.92±0.71) x 10~ 5 eV 2 ). For the 



le scale parameter £ of the distribution fl26|) to: 
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FIG. 6: The expected T2K sensitivity contour at 3 a CL, with two decoherence parameters con- 
tributing to the combined time evolution of gravitational MSW (stochastic metric fluctuations) 
damping with no energy dependence. 

completeness, we mention that, our best expected bound on the inverse-energy decoherence 
will imply, according to (I29p the bound on the life time t v jm v > 3 x 10 22 GeV -2 . 

The precise energy and length dependence of the damping factors is an essential step 
in order to determine the microscopic origin of the induced decoherence and disentangle 
genuine new physics effects from conventional effects, which as we have seen in section IHH 
may also contribute to decoherence-like damping. Some genuine quantum-gravity effects, 
such as MSW like effect induced by stochastic fluctuations of the space-time, are expected 
to increase in general with the energy of the probe, as a result of back reaction effected 
on space-time geometry, in contrast to ordinary-matter-induced 'fake' CPT violation and 



'decoherence- looking' effects, which decrease with the energy of the probe [49|]. At present, 
as one can see from the section IIII[ the sensitivity of the experiments is not sufficient to 
unambiguously determine the microscopic origin of the decoherence effects, but according to 
our estimations of the most plausible energy-length dependencies for the MSW like decoher- 
ence the sensitivity of CNGS and T2K will improve the current limits by at least two orders 
of magnitude and one would arrive at definite conclusions on this important issue. Thus 
phenomenological analyses like ours are of value and should be actively pursued when the 
data from OPERA and T2K will become available. When the present paper was finished we 
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FIG. 7: The expected T2K sensitivity contour at 3 a CL, with two decoherence parameters con- 
tributing to the combined time evolution of gravitational MSW (stochastic metric fluctuations) 
damping which are proportional to the neutrino energy. 



got aware on a similar analysis 59] performed for J-PARC experiments which agrees with 
our results concerning T2K and T2KK. 

In general, the characteristic energy dependencies of damping features are very interesting 
to search for physics beyond the standard model. In some cases, such damping signatures 
could be compensated by a shift of the neutrino oscillation parameters, which means that 
given such a damping effect, it is quite possible to obtain an erroneous determination of 
these parameters. However, if the damping effects are strong enough, then an establishment 
of effects beyond the standard neutrino oscillation scenario will be possible. Once such a 
damping effect is established, it will be very interesting to know from which non-standard 
mechanism it actually arises. Given this identification problem, we have found quite a low 
sensitivity for the models with inverse energy dependence in the exponent, which means 
that the damping effects of such kind are strongly correlated with the standard neutrino 
oscillation parameters, i.e., it is difficult to distinguish them from small adjustments in the 
oscillation parameters at CNGS and T2K. However, damping signatures similar to energy 
dependence or energy dependence squared can be very easily disentangled from the standard 
oscillations, but it is difficult to distinguish them from each other. Concerning different time 
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dependencies in the exponents including the combined signatures we have analysed, it can, 



in principle, be resolved if there are two baselines, as applied for example in [59J, and all 
the other parameters are known. Also, for a specific model, there may be relations among 
different 7's in ( 1441) and (|47[) like fits that actually imply much fewer independent parameters. 
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APPENDIX A: SIMMULATION OF NEUTRINO BEAMS 

A wide band neutrino beam is produced from the decay of mesons, mostly 7r's and K's. 
Mesons are created by the interaction of a proton beam into a needle shaped target, they are 
sign-selected and focused in the forward direction by two large acceptance magnetic coaxial 
lenses, conventionally called at CNGS (CERN) horn and reflector, and finally they are let 
to decay into an evacuated tunnel pointing toward the detector position. 

In case of positive charge selection, the beam content is mostly from the decay of 
7r + and K + . Small contaminations of I7 M (from the defocused 7r~ and K~) and v e (from 
three-body decay of K's and /x's) are present at the level of few percent. 



The neutrino fluxes for such a kind of beam are relatively easy to predict [53J once the 
secondary meson spectra are known, because the meson decay kinematics is well understood 
and the geometry of the decay tunnel is quite simple. 

Uncertainty in the estimation of the neutrino fluxes could arise because secondary mesons 
are selected over a wide momentum range and over a wide angular acceptance (~ 20 mrad). 

Re-interactions of secondary mesons in the target and downstream material contribute 
to reduce the neutrino fluxes and increase the uncertainty in the calculations (mainly for 
the wrong sign and wrong flavour contaminations). They are generally minimized using a 
target made of a number thin rods of low Z material interleaved with empty spaces (to let 
the secondary mesons exit the target without traversing too much material). In addition the 
amount of material downstream of the target (i.e. horn and reflector conductor thickness) 
is kept to the minimum. 
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The parameterization of the secondary meson production from protons onto a thin target, 
proposed in [53], is thus well suited to be used in neutrino beam simulations both because 
it extends its prediction over a wide range of longitudinal and transverse momenta and 
also because the small fraction of tertiary production from re-interactions in the target 
and downstream material can be accounted for with the approximations described in 531. 

n 

A comparison of the neutrino flux prediction based on the parameterization of [53J with 
some measured spectra is thus an effective estimator of the quality of the secondary mesons 
parameterization. For this purpose, the parameterization [53J has been coupled with a 
neutrino beam simulation program to be able to provide rapid and accurate predictions 
of neutrino spectra at any distance (i.e. short and long base line). The comparison has 
been performed both with already published data (CHARM II) and with predictions for the 
future CNGS long-baseline neutrino beam generated with GEANT and/or FLUKA based 
Monte Carlo programs. 



The resulting code 58] is a stand-alone application developed on the basis of parametriza- 
tion [53|] that allows to vary and optimize all elements and the geometry (in 3-D) of the beam 
line providing the results in terms of neutrino spectra and distributions at large distance 
with high statistics and in short time. 

The underlying idea is that in order to produce rapidly a neutrino spectrum at large 
distance over a small solid angle at CNGS beam, one has to force all the mesons to de- 
cay emitting a neutrino, and force all neutrinos to cross the detector volume. A weight is 
then assigned to each neutrino, proportional to the probability that this process actually 



happened. In practice this method is imp 
subsequent steps, as described in detail in 



emented by subdividing the simulation into four 



531 ] . These steps include: mesons production pro- 



cesses along target, meson tracking in the neutrino beam-line, neutrino production processes 
from mesons and neutrino production from muons. The weights associated with every step 
is described in details in [53]. 
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